Dynamics of a magnetic moment induced by a spin-polarized current 
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Effects of an incoming spin-polarized current on a magnetic moment are explored. We found 
that the spin torque occurs only when the incoming spin changes as a function of time inside of 
the magnetic film. This implies that some modifications are necessary in a phenomenological model 
where the coefficient of the spin torque term is a constant, and the coefficient is determined by 
dynamics instead of geometrical details. The precession of the magnetization reversal depends on 
the incoming energy of electrons in the spin-polarized current. If the incoming energy is smaller 
than the interaction energy, the magnetization does not precess while reversing its direction. We 
also found that the relaxation time associated with the reversal depends on the incoming energy. 
The coupling between an incoming spin and a magnetic moment can be estimated by measuring the 
relaxation time. 

PACS numbers: 75.70.Cn,72.25.Ba,75.60. Jk 
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I. INTRODUCTION 

Tremendous attention has been paid to the dynamics 
of magnetization in recent years because this problem is 
of fundamental importance in understanding magnetism 
and because the problem is of interest to technological 
applications in magnetic devicesi One of intriguing fea- 
tures of magnetization motion is spin transfer from a 
spin-polarized current to a magnetization of a ferromag- 
netic film, theoretically proposed by SlonczewskiS and 
Berger— , and later experimentally verified. 4,5 Since this 
spin transfer mechanism was first conceptualized, many 
studiesSi2i2ii£ have been performed on this phenomenon. 
However, the dynamics of a magnetic moment driven by 
a spin-polarized current has not been fully explored. 

In this paper we investigate the current-driven preces- 
sion and reversal of a magnetic moment. This is done 
quantum mechanically using a simple Hamiltonian, with- 
out introducing an external magnetic field. In this way 
one can easily distinguish contributions from the current 
from those induced by an externally applied field. To 
this end, we describe the motion of a magnetic moment 
in the lab frame where details of the magnetization re- 
versal are best illustrated. Since dynamics of a magnetic 
moment can be formally described in the local moment 
frame and such a description may also give some intuition 



J 



about the dynamics, we examine an interaction between 
a spin-polarized current and a magnetic moment in the 
local frame at the Hamiltonian level in section (II) . Then 
dynamics in the lab frame is illustrated in section (III). 
In this section one can see details of the dynamics such as 
under what conditions the motion of the magnetization 
can be non-precessional or the relaxation time associated 
with the reversal is a minimum. These phenomena have 
not been explored in the literature so far. Section (IV) 
is devoted to discussions about the adiabatic approxima- 
tion used to describe the motion of a magnetic moment, 
and we close with a summary. 

II. FORMALISM IN THE LOCAL MOMENT 
FRAME 

To describe effects of an incoming spin current on a 
magnetic moment M, as in Ref£ one can choose a frame 
(X'Y'Z'), where z' is parallel to M. Such a frame is 
called the local magnetic moment frame. Extensive work 
on ferromagnetism in the local moment frame has been 
done in Refiii An advantage of this frame is that it is 
trivial to diagonalize an interaction between an incoming 
spin s and a magnetic moment: — 2 Jff M ■ s, where Jh is 
the coupling. 

Let us start with a simple Hamiltonian relevant to the 
interaction: 



X 



H 



dx 



2m 



1> a {x) - 2J H M(x) ■ s + V{x)rb+{x)^ a {x) 



(1) 



where 4>a(x) creates an electron with a spin a at x, m 
is the electron mass and V(x) is an impurity potential. 
The magnitude of the magnetic moment is Mo, which re- 
mains unchanged. The electron spin can be represented 
as s l — (l/2)ip+a t ipi3, where a 1 is a Pauli matrix with 
i = x,y and z. We assume that the magnetic moment 



is determined by a localized electron \1/ so that the ki- 
netic part of the localized electron is not included in the 
Hamiltonian. 

Suppose a local magnetic moment M(x) points in the 
direction (6, <fi) at x as seen in Fig. 1. Then, a local rota- 
tion ( or coordinate transformation to the local moment 



2 



frame) is introduced: ip a (x) = U a p(x)xp{x) , where 



U(x) 



cos(6»/2)e"^/ 2 -sin(6»/2)e-^/ 2 
sin(0/2)e^/ 2 cos(6>/2)e 1 ^ 2 



In terms of x( x )i the Hamiltonian can be written as 



H = dx 



[xppo) ■ V (f/ Q7 x 7 ) - J ff x+C/+ Q (M • a afi )U^ X u + V{x)xtx a 



2m 

Since the interaction term in the Hamiltonian is diagonalized in this x( x ) basis, we obtain 



(2) 



(3) 



where 



H = H n 



dx 



A Q /3 • jap + A^pp a p 



(4) 



H 



ia/3 = 

AO) _ 

A ap — 



dx 



2m 



Xa ~ J H M X a O- al3 Xf3 + V(x)xt (x)Xa(x) 



— (VC7+ 7 ) • (VU 7 p) , and p a p = xtxp ■ 



r 



After diagonalizing the interaction, we have an extra 
term H ' = J dx A a/ 3 ■ j Q) g + A^lp a p in Eq. J3J instead 

of off-diagonal terms of the interaction in Eq. JJ). Using 
the explicit form of U(x), we can calculate vector po- 
tentials A^g and A a p. This was the route followed in 
Ref£, which led to a monopole-like term in the energy. 
Those authors attributed the spin torque term to this 
new vector potential, which is purely geometrical. 

Here we follow a different route, since we are interested 
in a simpler case, where the magnetization is not a func- 
tion of position. Thus, in our case of a single-domain 
ferromagnet, the extra term shown above will disappear 
because V£7 = 0. Instead, our spin torque will be present 
due to the dynamics of the coupled spin-moment system. 
In addition, we will not require an assumption regarding 
the magnitude of Jr in order to proceed, and we will uti- 
lize an impurity potential for convergence purposes which 
is otherwise irrelevant to the spin transfer as in Ref & 



III. DYNAMICS OF A MAGNETIC MOMENT 
IN THE LAB FRAME 

A disadvantage of the description in the local moment 
frame is that the precession of the magnetic moment can- 
not be seen; in other words, a precessional reversal of the 
magnetic moment cannot be distinguished from a plain 
reversal. Since our goal in this paper is to investigate 
the dynamics of the magnetic moment as mentioned in 



the introduction, we describe the motion of the magnetic 
moment in the lab frame. The geometry of our problem 
is shown in Fig. 1. We assume a single-domain ferro- 
magnet in the YZ plane for simplicity and consider the 
Hamiltonian Eq.Q. The incoming spin is along z and 
the direction of the magnetic moment is defined by 9(t) 
and 4>{t), which vary as functions of time t. 

The equation of motion for the magnetic moment 
M can be obtained quantum mechanically: dM/dt — 
i [H,M]. Since M i = (l/2)7 *+r* /3 * /3 , where * and t 1 
are the operator and a Pauli matrix for localized elec- 
trons, respectively, and 70 is the gyromagnetic ratio, the 
equation becomes 



dM 
~dT 



2j J H (M x s) 



(5) 



To analyze this equation we consider M as a classical 
vector and take s as its expectation value over the fer- 
romagnet. If we decompose s into a parallel sn and a 
perpendicular sj^ component to M, we know that only 
sj_ contributes to the equation. We can express sj_ using 
any unit vector. Let us choose, for the unit vector, the 
initial direction of the incoming spin §0 = z. Then 



s_l = s_l(M x s ) + s' ± M x (so x M) 



where s 



M-xy) and s' 

l-(So-M) 2 -L 



Iq-s-(J -M)(s-M) 



l-(so-A^) 2 

Eq. jnj, we can rewrite Eq. 10) as follows: 

dM 
~dT 



(6) 



Using 



-2 l0 J H s ± Mx(s xM)+2 l0 J H s' ± (Mxs ) . (7) 



3 



in/ 



're/ 




-> X 




-J H M o 



FIG. 1: Geometry of a quantum mechanical problem asso- 
ciated with the spin transfer. The incoming electron to the 
positive X axis are spin-polarized along z axis. The ferromag- 
net surface is at x — and parallel to YZ plane. The direction 
of the magnetic moment is defined by 9 and <f>, which are func- 
tions of time t. The ferromagnet is assumed to be sufficiently 
thick. 



As we can see in the above equation, the first term on 
the right hand side gives the spin torque while the second 
term causes a precession of the magnetic moment. We 
emphasize that the spin torque occurs only when s(t) 
changes as a function of time t. If s remains parallel to 
so, then s± vanishes and no spin torque takes place. In 
this instance, the effect of a spin is the same as that of 
an external magnetic field along z and the magnetic mo- 
ment precesses. In a phenomenological model^ the spin 
torque is represented by M x (so x M) with a proportional 
constant. However, a time dependence of sj_ is crucial as 
we emphasized. We also should stress that s± and s' ± 
are determined by dynamics, not geometrical details as 
in Ref-ifi 

To evaluate the expectation value of s, we need to solve 
the Schrodinger equation for the Hamiltonian Eq. QJ. 
Basically, the equation is one-dimensional because of 
translational symmetry in the YZ plane. We choose the 
direction of the polarized spin to be z. Then, an incom- 
ing wave function \ipi n ) with a momentum k or an energy 
e = k 2 /2m is \+)e lkx , where |+) is the spin- up state in 
the lab frame. We need to consider a normalization factor 
C for \ipi rl ). Since this wave function describes an elec- 
tron beam, |C| 2 is the number of electrons N e per unit 
length in one dimension. Intuitively, the more electrons 
are bombarded into the ferromagnet, the stronger is the 
effect of spin transfer. We thus expect the time scale 
for the reversal to scale inversely with N e (the more the 
number of electrons, the faster the moment responds). 
Similarly, the time scale will be proportional to the mag- 
nitude of the local spin, Si oca i (= Mo/70) (the larger the 
moment, the longer it will take to reverse it). 



An energy band and relations among k 2 /2m, 



FIG. 2 

fc^|j/2m, and JhMq 
k 2 /2m > J H M . 



In this figure, it is assumed that 



The reflected (\ip re )) and transmitted {\iptr)) wave 
functions are eigenstates |xj) and \xi) of the interaction 
2J H M • s = J H M ■ a ; namely, J H M • cr|x T ) = JifM |x T } 
and • cr\xi) = — JhMq\xi)- Therefore, 



IVve) = ^tIxtXxtI 



Ri\xi)(xi\ 



-ikx 



(8) 



while 



TtIxtMxtI+K^ + tJxxKxJ+K^ , (9) 



where fej = y/k + 2m JjjMo and fc[ = \Jk — 2mJijMo as 
depicted in Fig. 2. If the energy of the incoming electron 



fc 



iki becomes pure imaginary 
k, and its corresponding wave 



is less than JhMq, 
where k± = \j2m JhMq 
function decays exponentially; e~ Kl ' x . 

For x < 0, \tp(x < 0)) = \i/j in ) + \ip re ) and for x > 0, 
\ip{x > 0)) = \iptr)- The coefficients -R^q) and are 
determined by matching conditions of wave functions and 
their derivatives at x — 0: 



R 



k-k 



HI) 



HI) 



k + k 



and T- 



2k 



HI) 



HI) 



k + k 



HI) 



(10) 



Note that we take \4>in) — \+)e lkx in the above deriva- 
tions. This means that the number of electrons in 
the incoming beam N e is unity for simplicity; however, 
when we numerically solve the equation of motion for 
a magnetic moment, we can control this parameter. In 
the Hamiltonian Eq. we also have an impurity po- 
tential V(x). We shall introduce mean free paths Zf 
and li for each channel due to the impurity, and as in 
Ref£ they serve as convergence factors such as e x / l " i 
and e~ x l l v when we average the expectation of s using 
\ip(x > 0)) over the ferromagnet. We assume that the 
thickness of the ferromagnet (L) is much larger than the 
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mean free paths: L 3> ^m- One may wonder if the 
matching coefficients change when the convergence fac- 
tors are introduced. They do change as, for example, 
fc-f — > fcf + i/Z-j-; however, the conclusions we make later 
remain unchanged as we verified. 

Now we can calculate the expectation value of s within 



the ferromagnet; (s l ) — (l/2)(ip t r\o' l \'4>tr) with i = x ,y 
and z. The average values of the expectation values are 
evaluated as (s l ) = (1/2) J Q L dx{iptrW l \^tr)- After some 
straightforward algebra, we obtain for incoming energy 
greater than JrMq 



J 



(8*) 



^-Rc[a* 7 ] + ^Re[/3*<5] + Re 



(3*-f + 5*a 



(sv) = -^Im[ 7 *a]-^Im[r/3]-Im 



(l// T + l//i)-i(fc T -fci) 



^-(H 2 -|7| 2 ) + ^-(l/3| 2 



Re 



(1/Z T + - i(fc T ~ fcj 



(11) 
(12) 
(13) 



where a = (1/2)T T (1 + m z ), f3 = (1/2)71(1 - m x ), 7 = 
(\/2)T^(m x + im y ), and 5 = —{l/2)T^{m x + im y ). Here 
m (= M / '70 S \ocai) is the unit vector of the magnetic 
moment; namely, m z = cos(8) and m x +im y = sin(0)e 1 ^. 

In our treatment, the incoming energy e = k 2 /2m 
is a control parameter and JhMq is a scaling param- 
eter. Experimentally, e can be controlled by adjusting 
the applied voltage while JhMq is uncontrollable be- 
cause Jh is a microscopic parameter. If e = 77 Jh Mo, 
then k 2 /2m = (t] + 1)J h M and kf/2m = (77 - 1)J H M . 
Defining k^/2m = JhMq, fcj and fc| can be written as 
fcf = V 7 ? + 1 ^0 and fcj. 



-^77 — 1 ho- Since the current 
density is in energy units in ID (H = 1), using j = k /m 
with one electron per unit length we can define a di- 
mensionless time r = j t, which will be used in the nu- 
merical calculations. When jj < 1, fej, = — r] ko 



as mentioned earlier 
k 



In this case (s) changes to reflect 



i\f\ — 77 kg. We do not present equations for 77 < 1 
here because the derivation is parallel to the above case 
and expressions are similar with those for 77 > 1 . Since we 
attribute the impurity potential to the mean free paths, 
it is natural to assume Zj = l[ = I. We also introduce 
a parameter a — lk - In the numerical calculations, we 
vary a from 0.5 to 2. Qualitative behaviors of m are not 
sensitive to the value of a. 

A dimensionless equation of motion for the magnetic 
moment is 



N e /2 

— = (m x h) 



dm 



S, 



ocal 



where 



(14) 



hi = \\T,\ 



; (1 



m z )m i --\T l \ 2 (l 



m z )mi 



4/a 2 + (^PTl-Vv^) 2 



(15) 



with (i = x, y, and z) 



A x 
B x 
A v 



(Re [T T Tf] m x m z + Im [T T Tf] m v ) 



(Re [T|T*] m y - Im [T T Tj 



I J m x m z ) 



(Re [7\T*] m v m z - Im [T T T X *] m x ) 



B v = ^ (Rc [T T Tf] m x + Im [T T Tf] m y m z ) 



A z = -Re [T-fT*] (ml+m 2 y ) 



B z = -Im[T T Tf] (m 2 x + m 2 y ) . 

Clearly the factor N e /2Si oca i could be absorbed into the 
time (already dimensionless). Since its effect is obvious, 
we set N e /2Si OC ai — 4 for all our results. 

We choose various values of 77 between 0.25 and 4, and 
show m,(r) vs. r and a locus of m in the (m X} m y ,m x ) 
coordinate. For an initial condition of m we choose 
9q = 7r/1.01 and <fio = ir/4 to see the magnetic moment 
reversal. Because of a rotational symmetry, the initial 
value of (f> is not important. It is obvious that if 9 = 
or 7r, the spin polarized current has no effect on m. In 
Fig. 3(a), we show the locus (dotted curve) of m for 77 = 2 




FIG. 3: Precessional reversal of the magnetic moment for 
r\ — 2 and a — 1. Fig. 3(a) shows the locus (dotted curve) of 
m and Fig. 3(b) is for mi("r) vs. r. The initial direction of m 
is given by do — 7T/1.01 and (po — ir/4. Thin circles define a 
uni-sphere. 



FIG. 4: Plain reversal of the magnetic moment for r/ — 0.5 
and a = 1. Fig. 4(a) shows the locus (dotted curve) of m and 
Fig. 4(b) is for mj(r) vs. r. The initial direction of m is the 
same as in Fig. 3. Note that there are no oscillations in m x 
and m y . Thin circles define a uni-sphere. 



and o = l, and plot mj(r) vs. r in Fig. 3(b). Thin cir- 
cles define a uni-sphere. Oscillations in m x and m y imply 
precession of m. For r/ = 2, m shows a precessional rever- 
sal. On the other hand, for r\ — 0.5 it has a plain reversal 
without precession as we can see in Fig. 4(a) and (b). In 
this instance, m x and m y do not show oscillations. The 
precessional reversal takes place only when rj > 1. This 
remains true for a = 0.5 or 2. We plot these results in 
Fig. 5(a) and 5(b) for r\ = 0.25 and 4. 

One can define the relaxation time tq of the reversal 
as an elapsed time during the reversal between 9 ~ ir 
and 6 ~ 0. When m z ~ 1, we can parametrize 
In [1 — m z (r)] = C\ — c 2 t/t , where C\ w 8.9 and C2 ~ 
13.6. We found these values are independent of r\ and 



a. For given 77 and a, we can determine r by comparing 
numerical results with c\ — C2t/tq. For example, To — 7.9 
for 77 = 0.9 and o = 1. In general, the smaller a (or I) is, 
the longer r is for a given 77. This can be understood be- 
cause the wave function \iptr) decays faster if I is shorter 
so that the spin transfer is relatively less effective and, 
thus, it takes a longer time to reverse m. In Fig. 6, we 
plot m 2 vs. t for i] = 4 (main frame) and for rj = 0.25 
(inset) with a = 0.5 (solid) 1 (dashed), and 2 (dotted 
curve). In this figure, we can see the relation between 
a and tq mentioned above. For r] > 1, as a increases, 
weak precession occurs because r decreases as seen in 
the main frame of Fig. 6; in other words, m does not 
have enough time to precess strongly. We can also see 




FIG. 5: Locus of m for r\ — 0.5 (thick dotted curve)and r\ — 4 
(dotted curve). In Fig. 5(a), a = 0.5 while a = 2 in Fig. 5(b). 
Regardless of a, no precession occurs when r\ = 0.25. Thin 
circles define a uni-sphere. 



such a behavior in Fig. 5 comparing a — 0.5 and a = 2 
for r] = 4. 

We plot To vs. -q in Fig. 7 for a given a. The re- 
laxation time is evaluated using the parameterization: 
In [1 — m z (r)] = c% — citJtq. In the main frame, a = 1 
while in the inset a = 2. At r = to, m z (ro) ~ 0.99 for 
all plots. Interestingly, tq is minimum at r\ ~ 1. There- 
fore it is possible to estimate the microscopic coupling 
parameter Jh between an incoming spin and a magnetic 
moment by measuring tq{t)), because r has a minimum 
for a given mean free path. 



IV. DISCUSSION AND SUMMARY 

In this section we would like to discuss the adiabatic 
approximation, which we tacitly used to study the dy- 
namics of a magnetic moment. First we summarize 



FIG. 6: m z as a function of r. In the main frame, 77 = 4 while 
in the inset r\ — 0.25 with a = 0.5 (solid) 1 (dashed), and 2 
(dotted curve). 




4 1 1 1 1 I 

o.s 0.9 1 1.1 1.2 " 

1 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 

0.5 1 1.5 2 2.5 

FIG. 7: The relaxation time to vs. r\. In the main frame, 
a — 1 while in the inset a — 2. to has a minimum value at 

7,-1. 



the procedure we followed. We calculated (s) using 
\t()(x)); namely, (s) = (l/2)(tp(x)\a\tp(x)) for x > to 
solve <IM./dt — 2jqJh (M x (s)). Here we mention that 
\ip(x > 0)) is obtained by considering the Hamiltonian at 
a given time t following Ref£ Since the incoming wave 
function \ipi n ) ~ |+) is not an eigenstate of the Hamilto- 
nian for x > 0, we have a linear combination of |+) and 
|— ) for \ip(x > 0)) and \ip(x < 0)). The matching condi- 
tions of wave functions at x = allow us to express the 
coefficients of the combination for \i/j(x > 0)) in terms of 
M(i) (see Eqs. © and JXDJl). Now (s) is a function of 
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M(t), and the time dependence of (s) is given exclusively 
by M(t). This means that the time evolution of the wave 
function for x > is not fully taken into account. In ad- 
dition to the equation for dts/L/dt, one can derive the time 
derivative of the spin operator using ds/dt = i [H, s] : 

ds 

— + V • J = 2 J H (s x M) , (16) 

where J is the spin-current tensor. It is obvious that 
when we calculate an expectation value of s in Ea.l |16(l 
we need to use \ip(x,t)); (s) t = (l/2)(i/)(x,t)\cr\ilj(x,t)), 
where \ip(x,t)) is obtained from i-^\ip(x, i)) = H\tp(x,t)}. 
Rigorously speaking, one has to solve the two coupled 
equations for M and s using \ip(x,t)) to calculate the ex- 
pectation value of s and V ■ J . However, if we compare 
Eq. © or (T2J with Eq. JTEJl, we see that Eq. has 
a factor l/Si oca i while Eq. (|16f) does not. This means 
that if we treat the magnetic moment semiclassically, i.e. 
Siocai 3> 1, then the time scale of Eq. Q14f) is much longer 
than that of Eq. (|16|l . Therefore, the adiabatic approxi- 
mation is applicable to our analysis. 

In summary, we have studied the effect of an incoming 
spin-polarized current on a local magnetic moment in a 
magnetic thin film. We found that the spin torque occurs 
only when the incoming spin changes as a function of time 
inside of the magnetic film. If the incoming spin remains 



parallel to its initial direction, no spin torque takes place. 
This implies that some modifications are necessary in a 
phenomenological model where the coefficient of the spin 
torque term is a constant. Moreover, the coefficient is 
determined by dynamics instead of geometrical details. 
The magnetization reversal can be precessional as well 
as non-precessional depending on the incoming energy of 
electrons in the spin-polarized current. If the incoming 
energy is greater than the interaction energy (JjjMo), 
the magnetization precesses while reversing its direction. 
For the incoming energy smaller than JhMq, the magne- 
tization reversal is non-precessional. We also found that 
the relaxation time associated with the reversal depends 
on the incoming energy for a given mean free path. Our 
numerical calculations imply the coupling between an in- 
coming spin and a magnetic moment Jh can be estimated 
by measuring the relaxation time. 
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